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Abstract 

We are concerned with the acoustic scattering problem, at a frequency k, by many small obstacles of 
arbitrary shapes with impedance boundary condition. These scatterers are assumed to be included in a 
bounded domain Q, in which is embedded in an acoustic background characterized by an eventually 
locally varying index of refraction. The collection of the scatterers Dm, rn = is modeled by 

four parameters: their number M, their maximum radius a, their minimum distance d and the surface 
impedances Am, m — We consider the parameters M,d and Am’s having the following scaling 

properties: M := M{a) — 0{a~‘’), d := d{a) ~ a* and Am := Xm{a) = Am.ou”^, as a —>■ 0, with non 
negative constants s,t and fd and complex numbers Am,o’s with eventually negative imaginary parts. 

We derive the asymptotic expansion of the farfields with explicit error estimate in terms of a, as a —>■ 0. 
The dominant term is the Foldy-Lax field corresponding to the scattering by the point-like scatterers 
located at the centers Zm's of the scatterers Dm’s with \m\dDm\ as the related scattering coefficients. 
This asymptotic expansion is justified under the foliowing conditions 

a < ao, |K(Am,o)| > A-, |Am,o| < A+, /3<1, 0<s<2-^, 

and the error of the approximation is C as a —>■ 0, where the positive constants ao, A-, A+ 

and C depend only on the a priori uniform bounds of the Lipschitz characters of the obstacles Dm’s 
and the ones of M{a)a‘ and . We do not assume the periodicity in distributing the small scatterers. 
In addition, the scatterers can be arbitrary close since t can be arbitrary large, i.e. we can handle the 
mesoscale regime. Finally, for spherical scatterers, we can also allow the limit case /3 = 1 with a slightly 
better error of the approximation. 
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1 Introduction and statement of the results 

Let Bi, B2, ■ ■ ■, Bm be M open, bounded and simply connected sets in with Lipschitz boundaries con¬ 
taining the origin. We assume that the Lipschitz constants of Bj, j = are uniformly bounded. 

We set Dm ■= eBm + Zm to be the small bodies characterized by the parameter e > 0 and the locations 
Zm € m = 1,. .., M. Let W be a solution of the Helmholtz equation (A -h k^)U'' = 0 in We denote 
by C/® the acoustic field scattered by the M small bodies Dm C due to the incident field C/* (mainly the 
plane incident waves U'^{x,9) := with the incident direction 9 G S^, where being the unit sphere), 
with impedance boundary conditions. Hence the total field C/‘ := C/* -I- H® satisfies the following exterior 
impedance problem of the acoustic waves 

Dr^ , ( 1 . 1 ) 
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— inU'^ = o 



0, 1 < TO < M, 


|a;| —)■ 00 , (S.R.C) 


(1.3) 


where Vm is the outward unit normal vector of dDm, k > 0 is the wave number and S.R.C stands for the 
Sommerfeld radiation condition. The scattering problem ( 1.1|1.3 1 is well posed in the Holder or Sobolev 
spaces, see 18 19 27 ^ for instance, in the case > 0. We will see that this condition can be relaxed 
to allow ’XsXm to be negative under some conditions. Applying Green’s formula to C/®, we can show that the 
scattered field U^{x,9) has the following asymptotic expansion: 


u%x,e) 


giK.\x\ 

47r|a;| 


c/~(i,0) + o(|x|-2). 


|a:| —)■ 00, 


(1.4) 


with a; := 1 ^, where the function U°°{x, 9) for {x, 0) € x is called the far-held pattern. We recall that 
the fundamental solution, $^(3^,2/), of the Helmholtz equation in with the hxed wave number k is given 
by := for all X, y e 

Definition 1.1. We define 

1. a ■.= max diam(Dm) \ = e max diam(Bjn)], 

l<m<M l<m<M ^ 

2. d := min dmj, where dmj '■= dist{Dm, Dj). 

3- Kmax as the upper bound of the used wave numbers, i.e. k G [0, Kmax]- 


I 


The distribution of the scatterers is modeled as follows: 

the number M := M(a) := 0(a“®) < MmaxO,~^ with a given positive constant M^ax 


5. the minimum distance d 

dmin and dmax- 


d{a) « a*, i.e. < d{a) < djnaxoX, with given positive constants 


6. the surface impedance Am := Xm,oa ^, where Xm,o 7^ 0 and might be a complex number. 

Here the real numbers s, t and j3 are assumed to be non negative. 

We call the upper bounds of the Lipschitz character of Hm’s, Mmax, dmin, dmax and Kmax the set of the 
a priori bounds. 

The goal of our work is to derive an asymptotic expansion of the scattered held by the collection of the 
small scatterers Dm, rn = 1, •.., M, taking into account these parameters. This is the object of the following 
theorem. 


Theorem 1.2. There exist positive constants ao, A_ and A+ depending only on the set of the a priori bounds 
such that if 

a < oo, |Am.o| < A+, |3fi(Am.o)| > A_, /3<1, s<2-/3, (1.5) 

then the far-field pattern U°°{x,9) has the following asymptotic expansion 

M 

U^{x, 9)=J2 + o (a3-®-2/3) , 

m—1 


( 1 . 6 ) 














uniformly in x and 9 in S^. The constant appearing in the estimate 0{.) depends only on the set of th^ 
a priori bounds, A_ and A+. The coefficients Qm, Tn = are the solutions of the following linear 

algebraic system 
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Qm “t" 'y ^ ^j)Qj — ^mU i 


(1.7) 


1=1 

j^m 


for m = 1, M, with Cm ■= -Xm\dDm\- 

The algebraic system o* is invertible under the conditins: 


|3?(Am.o)| > A_, s<2-/3. 


( 1 . 8 ) 


Let us now assume that the background is not homogeneous but modeled by a locally variable index of 
refraction n, i.e. there exists a bounded set ft such that n{x) = 1, a; G \ fl and bounded inside fl, i.e. 
|n(a;)| < Umax, a; G n. In this case we model our scattering problem as follows. The total field := W + Uf, 
satisfies the following exterior impedance problem of the acoustic waves 


(A + K"^n^(a;))?7„ = 0 in IR'^\ U D 
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Again, the scattering problem ( 1.9|1.11|) is well posed in the Holder or Sobolev spaces, see 18 19 2^ in the 
case 9Am > 0. As we said for ( 1.9||1.11 ), this last condition can be relaxed to allow QXm to be negative. 
Applying Green’s formula to Uf, we can show that the scattered held Uf{x, 9) has the following asymptotic 
expansion: 

„iK.\x\ 

K{x,9) = ^^U^{x,9) + 0{\x\-\ |a:|^oo, (1.12) 

where the function Uf°{x,9) for {x,9) G x is the corresponding far-held pattern. As a corollary of 
Theorem 1 1.2[ we derive the following result: 


Corollary 1.3. There exist positive constants ag, A_, A+ depending only on the set of the a priori hounds 
and on Umax such that if 


a < oo, lAm.ol < A+, |3fJ(Am,o)| > A_, /3<1, s<2-/3, -<t 


then the far-field pattern Uf°{x,9) has the following asymptotic expansion Q 


M 

Uf^{x,9) = Vf^{x,9) + 5] Vf{zm-x) 

m—1 


+ O (a3-^-2/5) 


(1.13) 


(1.14) 


uniformly in x and 9 in S^. The constant appearing in the estimate 0{.) depends only on the set of the a 
priori bounds, A_, A+ and Umax- The quantity *) is the total field, evaluated at the point Zm in 

the direction —x, corresponding to the scattering problem 


(A -I- K^n^{x))Vf = 0 


(1.15) 

^ The following remarks apply for both Theorem El and Corollary o_ If s < 2 — /3, then the bounds A_ and A_|_ can be 
arbitrary. If KA^n < 0 then, for the invertibility of the systems in and | |l.l'7| >, the condition s < 2 — /3 on the number of the 
small bodies can be replaced by the condition t < 2 — ^ on the minimum distance between them, see Remark |6.1| 

^ By the mixed reciprocity relations, we have —^) = G^ {x,Zm) wh ere G^{x, Zm) is the farfield created by a point 

source GH^{x,z^n), located at the point Zm-, in the scattering model 1 1.15|1.16 i. 
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- inV^ = o 



The coejjicients Qm, m = 1, are the solutions of the following linear algebraic system 


M 

E 

i=i 

j^m 


CmGK,{z 


= -c^v\z^,e), 


(1.17) 


for m = 1 ,..., M, with 

Cm ■■= -Xm\dDm\. 

Here Gf^{x,z) is the outgoing Green’s function corresponding to the scattering problem (1. 
The algebraic system (1.171 is invertible under the conditions: 


(1.18) 


|3?(A^,o)| > A_, s<2-(3. 


(1.19) 


Before discussing these results comparing them to some of the literature, we add the following remark 
on the particular case when the scatterers are spherical. 


Rema rk 1.4. If the scatterers have spherical shapes, then we c an slightly improve the error estimates in 
{TI 4 ) by slightly changing the form of the algebraic system (1.17). Precisely we have O instead of 

7-s-2/3^_ For these spherical shapes, we can also handle the case (3=1. The result reads as follows. 


O 


There exist positive constants qq, A_, A+ depending only on the set of the a priori bounds and on Umax such 
that if 

a<ao, |A™,o| < A+, |5ft(A,„,o)| > A_, /3<1, s<2-/3, (1.20) 

then the the far-field pattern Uff{x,6) has the following asymptotic expansion 


M 

U^{x,e) = V(^{x,e) + ^ V:,{Zm,-mra+0{a^-^-^) , ( 1 . 21 ) 

m—1 

uniformly in x and 9 in S^. The coefficients Qm, rn = 1,...,M, are the solutions of the following linear 
algebraic system 
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Y. CmG^{: 


■m.'i 


= -CmV\zm,e), 


( 1 . 22 ) 


i=i 

j^m 


for m = 1,..., M, where now we have a slight change in the constant Cm compared to {1.18), i.e. 

^ \m\dDm\ 

■“ -1 4- X T ■ 

A- \ 

where Im ■= Jgjj ^o{sm,t)dsm, t G dDm, is a constant if Dm is a ball, rl 
This algebraic system is invertible under the same conditions as ^1.1 TV. 


(1.23) 


These results say that the dominant term in the approximation (1.6), and similarly the one in (1.14), is 
the Foldy-Lax field corresponding to the scattering by the point-like scatterers located at the ’centers’ Zm’s 


® The property that Im is a constant only for ba lls i s known as Gruber’s conjecture and it is solved for the 2D for Lipsc hitz 
regular domains and in 3D for convex domains, see [23| . This is the main restriction why we cannot obtain in Theorem |1.2| and 
Corollary 1 1. 3| the same results as in Remark |1.4| for general shapes. Actually, if Dm is a ball then the value of Im is nothing 
but the radius of Dm (= e radius(Bm)). 














of the scatterers D^s with Xm\dDm\ as the related scattering coefficients. This Foldy-Lax field describes the 
field that results in the multiple scattering between the different scatterers Dm’s, see 


22 and the references 


therein for more information on this issue. The accuracy of the approximation of the scattered field by the 
Foldy-Lax field depends of course on the error estimates. In its generality, this issue is still largely open but 



all the involved parameters: the maximum diameter a, the number M, the minimum distance d and the 
surface impedance A^’s including the regime defined by M = M{a) := a“®, d := d{a) := a‘, Xm '■= Xm(a) '■= 
Am.oa”^. We characterized the set of parameters s,t and /3 where the approximation makes sense and we 
provided the approximation with explicit error estimates in terms of these parameters. These approximation 
formulas can be used for different purposes: 

1. First, they can be used for imaging where the small anomalies are modeled by the small bodies. There 
is a large literature on the mathematical imaging of small anomalies, see for instance and the 
references therein. Compared to those results, we can allow the small bodies to be dense and very 
close, since their number M can be very large and the minimum distance d can be as small as we 
want. This last property means that we deal with the mesoscale regime, where ^ ^ 1 

and minj^m dist{z 
with 


24 25 


4.J, , 


mindist{zj ,Zm) 

is the minimum distance between the centers of the small scatterers, compare 


2. Second, these approximations can be used for the design of new indices of refraction. When the small 
scatterers are periodically distributed, this can be justified by homogenization, see 


l3l[i7popT . Using 


the approximations in (1.6) and (1.14), we do not need such periodicity in the distribution of the small 
scatterers. This observation has been already made in [31] with quite formal computations. Let us 
observe that in our approximations in Theorem |1.2| and in Corollary |1.3[ we allow the surface impedance 
to have negative imaginary parts. With this freedom of taking the surface impedance, we can generate 
a large class of indices of refraction with possible applications to the acoustic metamaterials. Details 
on this issue are reported in and [^. Let us emphasize that we provide those results with explicit 
error estimates. This gives more credits to the feasibility of the design process. 

3. Third, these approximations can be used to extract the values of the index of refraction n{x),x G D, 
from the farfields corresponding to few incident directions. The idea of adding small inclusions to 
deform the medium and then extract the field inside the support of the medium is already described 
and used in the literature, see for instance [^. The extraction of the index of refraction (or other 
coefficients) from the internal measurements, related also to the hybrid imaging methods, has recently 
attracted the attention of many authors, see for instance [7- 10 and the references therein. One of the 
main difficulties is to handle the possible zeros of the total field. What we propose is to deform the 
medium using multiple (and close) inclusions instead of only single ones. In this case, what we derive 
from the asymptotic expansions are the internal values of the Green’s function and not only the total 
fields. Finally, the values of the index of refraction can be extracted from the singularities of these 
Green’s function. Hence, we avoid the problems coming from the zeros of the internal fields. These 
arguments are reported in [^. 


To finish this introduction, we discuss the two ’extreme’ cases given by the Neumann and the Dirichlet 
boundary conditions. In the former case, Am’s are all zero. In the results stated above, we see that Am’s 
are not allowed to be simultaneously zero. Otherwise all the coefficients Cm, m = 1 ,..., M, will be zero and 
hence all the constants Qm, 
where Am = 0 


m = vanish and the expansion (1.6) will make no sense. The case 

m = 1, ...,M, is quite tedious but interesting. 


Technically, we see in (3.73), for instance, that in this case we need to go to the higher order in 


the expansion. Doing that requires quite tedious computations remembering that we are taking into 
account all the parameters describing the scatterers (M, a and d). 


































• This particular case is interesting because in this case the dominant coefficients are defined by matric& 
and not vectors (as the vector (Ci, (72, Cm) in the case where Am ^ 0). But this is not a surprise since 
the dominant term of the expansion of the far-fields are the Foldy-Lax fields modeling the interaction of 
the multiple point-like scatterers (given here by the ’centers’ of the scatterers). For hard scatterers, we 
talk about anisotropic interaction between the scatterers, see the book [22| for more information about 
this issue, contrary to the impedance case where the interaction is isotropic. Due to this ’anisotropic’ 
character of the dominant term, the equivalent medium (when we distribute a cluster of small scatterers 
with Neumann boundary conditions) is characterized by a divergence form Helmholtz model where 
the coefficient appearing in the higher order derivative is a matrix (defined by the (scaled) matrix- 
coefficients appearing in the dominant term of the expansions). 


These arguments need to be mathematically justified and quantified. The approximation in the Dirichlet 
case is discussed in 15 where it is justified under the condition that y/M — is bounded, by a constant 
depending only some a priori bounds, which means, in the scales we use here, that 0 < s < 2 — 2t, or t < 

Since the small obstacles are distributed in a bounded domain H, then we have the natural condition on 
their number M — 0{d~^), i.e. | < t. In the present work, this last condition is the only one we (naturally) 
impose on t. Hence, compared to 15 , we allow here the small scatterers to be as close as we want since t 


can be as large as we want, i.e. we cover the mesoscale regime. However, we believe that the conditions used 


15 can be improved to be comparable to the ones we impose here. 


The rest of the paper is devoted to prove Theorem 1.2 and Corollary 1.3 In section 2, we describe briefly 

The detailed proof of Theorem |1.2| is done in section 3 while 


the main steps of the proof of Theorem 1.2 


the one of Corollary |1.3| is given in section 4. The justification of Remark |1.4| is discussed in section 5. In 
the appendix, we deal with invertibility of the algebraic system (1.7) (and similarly for ( 1.17|)). 


2 A brief description of the proof of Theorem 1.2 


Let us here describe very briefly the main steps of the proof of Theorem 1.2 Firts of all, the scattering 
problem has a unique solution and it can be represented via single layer potentials 


^ f ( 

{x) = U\x) +/ ^f,{x,s)ajn{s)ds, X 

m=l -^ 90 ^ A 


M - 

U Dr, 

m—l 


( 2 . 1 ) 


where a := (cti, ..., ctm) satisfies the corresponding system of integral equations. We show that this system 

M 

of integral equation is invertible in the space L‘^{dDm) and that the solution of the scattering problem 

m=l 

is unique under some natural conditions on the eventual negative imaginary parts of the surface impedance. 
We divide the rest of the analysis into few steps: 

1. We have the following a priori estimate of the densities Cm’s. If s < 2 — /3 then \\(Jm\\L‘^{dD,r^) — Ca}~^ 
where C depends only on the Lipschitz characters of Bm, rn = 1,..., M. The delicate work here is to 
derive the precise scaling of the corresponding boundary integral operators in the appropriate boundary 
Sobolev spaces (taking into account the three parameters s,t and /3). 


2. From the representation U‘^{x) = foo s)am{s)ds, for x G 

using the above a priori estimates on Cm’s, that 
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M - 

U Dm ), we deduce, 

m=l 


C/“(x) = ^ 
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= E 


e-“"'V„(s)ds 


'dD„ 


-IKX- Z, 


'dD„ 


'am{s)ds + 


'dD„ 




m—l 
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= E 


-IKX- 2, 


+ 0 {KMa' 


3-/5 
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( 2 . 2 ) 


where := (Tm{s)ds. 

3. To estimate the terms Qm, we use the boundary conditions. For Sm € dDm, using the impedance 


boundary condition (1.2), we have 
dU* 


0 = 


dVr 


H“ ^mJJ (^m) — 


i(Sm) 


f ^ 


(Sm,s)crr„(s)ds + V' / {Sm, s)aj {s)ds 

JdDj OVm 


J = 
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dW 


+A 

m EL 


j=l-^OD 


(2.3) 


Integrating the above on dDm, we obtain 

9$. 


f ^m{,^7n)dSjyi + [ ( [ i^)ds+y:[ f ) crj(s)ds 

JdDrr, JdDrr^ \jdDrn -laDj \jdDrr, J 
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+Am 

It can be rewritten as 


j=l-’9D, \JdD^ 


^)dSm I ^m(^)ds — 


r dw 

laD,,^ dvm 


{Sjn)dSyYi I XmU {^Sm)dSri 


ldD„ 


2 Qm + 


9^0 

^Drr^ \JdDrn 


('5m;5)dSy7j j (7‘YYi{^s')dS 
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dDi KJdDr^ 


('^m; si^dSjYi (jj{s)ds 


= -.A 


= :B 


+ Aj 


, / / d>„(sm,s)dsm ) crm(s)ds+Am ^ / / $„(sm,2j)dsm ) o-j(s)(is 

\JaDm _/_^ \JaD^ J 


=-.c 


= :D 
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I {^^m)dSjYi / XmU {^Sr[ri)dSjYi -\- A -\- XjyiD , 


with 


V^aD 
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j^m 




an, \jan. 




(2.4) 


dSm^ am{s)ds 

(2.5) 

s) - ^f,{Sjn,Zj)]dSmj (Tj{s)ds. 

(2.6) 


Using the a priori estimate of (im’s, the singularities of the fundamental solutions d>„ and the har- 
monicity of $ 0 j we derive the estimates: 

A=-^Qm, Al' = 0(s:V-^), i? = o(^2«2^[A + 7]j, C = O (a^-^) , 


(2.7) 























( 2 . 8 ) 
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D=Y,^n{zm,z,)Qj\dD^\+0[2 
j^m 


,5-/: 


d2° 


7k + 


6k + 13 


and 


D' = 0 [2 


<P°‘ 


Ik + 


6k + 13 


(2.9) 


Here the parameter a is introduced to count the number of small scatterers surrounding a given and 
fixed one, see Fig[T] and the discussion before. From ( |2.4[ ), we obtain the approximation below 
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Qm ^ ^ ^j 

j^m 

= + o ((|A™| + K)Ka^) + A„0 (a^"^) 

+0 {k^o^-^) + O f2K2-^[— + 7]^ + A^O [2 


,5-/: 


d2c 


above system as 


Qn 
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IKsO-. 


^ F) / 3~/5 \ 

-'^Cj^^(zm,Zj)-^ +0 \ a^ ^ ■ 

j^m J \ / 


7k + 

6k + 13 

) • (2.10) 






we can rewrite the 

3-/3 \ 

tP°‘ ) 


(2.11) 


4. Let now the vector Y := (Yi, Y 2 ,..., Ym) be the solution of the Foldy-Lax algebraic system 


M 


Ym = ^ z,)Y,, m = l ,..., M. 


( 2 . 12 ) 


1=1 

j^m 


We show that the algebraic system (2.12) is invertible under general condition on s and /3 and derive 
an error estimate. Based on this error estimate, we deduce from (2.11) and (2.12) that 


^ / / 3~/5 \ \ 

Y^lQm-CmYml = O (m a^-^ ^ V 

m—1 ^ ^ ^ 


(2.13) 


5. The proof ends by plugging (2.13) and (2.12) in (2.2) and setting Qm := CmYm, m = 1 ,..., M. 


3 The detailed proof of Theorem |1.2 

3.1 The representation via layer potential 

We start with the following proposition on the solution of the problem ( 1.1|1.3 ) via the layer potential 
representation. 


























Proposition 3.1. Assume that the negative part of the imaginary part of Xm are small enough. In additioi^, 
suppose that is not a Dirichlet^ eigenvalue of the Laplacian in Dm, for m = Then [3/or 

m = 1,2,..., M, there exists am G L^(dDm) such that the problem (1.1 1.3) has one and a unique solution 
and it is of the form 


U*(x) = U\x) +'V' / s)o-„^(s)ds, a; G ( U Dm ] , (3.1) 


Proof of Proposition 3.1 We look for the solution of the problem ( 1.1||1.3 ) of the form (3.1), then from 
the impedance boundary condition (1.2), we obtain 
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'dDi dvj{sj) 
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i{s)ds+ ^ 


d^K.{sj,s) ^ 


i{s)ds 


m^j 
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+^3Y] I ‘^K{sj,s)amis)ds = - \ - XW^isj), ^Sj G dDj, j 

m=iJdD^ 0^j[Sj) 


One can write it in a compact form as (—5I + DL* + DK* + X{L + K))a = —{di, + with di, := 

idumj)^.,=i, A := {Xmj)f^,,=„ DL* := {DL*mX.i^i^ DK* := {DK*mA^.^„ L := and K := 


'"mj — )m,j 

iKmj)m,j=i^ where 


I, Identity operator m = j 


0, zero operator 


else 


DL*m, = 


VZ ,, m = j 


mj 

0 


else 


DK* = 


D*mj rn^j 
0 else ’ 


Q _ J m = j 

1 0 else 


Xmj — 


Am m = j 
0 else 


Lmj — 


Smj 'OV — j 
0 else 




= U^{si, ..., sm) '■= {U'^isi),..., W{sm)) and cr = tT(si,..., sm) ~ (cti(si), ..., aMism))^ ■ Here, 
for the indices m and j fixed, Smj is the integral operator acting as 


5, 


mj 


[aj){t)~ / $„(t,s)crj(s)ds, t&dDm, 

JdDi 


(3.5) 


and Vmi is the adjoint of the integral operator defined by. 


Dmj{<Jj){t) ■= f '"^'"'^y\^ orj{s)ds, tedD„ 
Jdl 


d^^,{t,s) 

IdDi dVm{s) 


(3.6) 


Then the operator Vmm • L'^idDm) —t L‘^{dDm) is adjoint of the double layer operator Vmm '■ 
L‘^{dDm) L'^{dDm), defined by, 

Dm7n{<ym){t) ■■= f (s)ds, t G dDm, (3.7) 

JdD^ 


^This last condition is satisfied for every k such that k < Mimax and a < 


^ji/ 2 , 1 - Here ji/ 2,1 positive zero 


of the Bessel function Ji/ 2 - 

^ The result of this proposition is valid regardless of the smallness of the obstacles Dm’s nor the conditions on M, d and 
Am’s. The emphasize is on the possibility to deal with surface impedance eventually having negative imaginary parts. 
















and the operator —\l + T^mm ■ L'^{9Dm) L‘^{dDm) is isomorphic and hence Fredholm with zero indelj^^ 
For m j, : L'^{dDj) L'^{dDjn) is compact, see 28 Theorem 4.1]p| 


Also notice that Smj ■ L‘^{dDj) —>• L'^{dDm) is compact. So, (—+ DL* + DK* + \{L + K)) : 

M M 

n L'^idDjn) n L^{dDjn) is Fredholm with zero index. We induce the product of spaces by the 

m—l m—l 

maximum of the norms of the space. To show that {—^1 + DL* + DK* + X{L + K)) is invertible it is enough 

_L nr* j- nrr* 

M 


to show that it is injective, i.e. (—5I + DL* + DK* + \{L + K))a = 0 implies cr = 0. We write 


U{x) = Y, 


-1 ^dDr, 


m—l 


s)CTm(s)ds, in K'^\ ( U D^ 

' m=l 


and 


M 


U{x) = Y 


M 


'dDrr 


^^{x,s)u^{s)ds, in U D^- 


m—l 


Then U satisfies AU-\-k‘^U = 0 for x € 


M - 


M 


U Drn , with S.R.C and +A^n)?7(x) = 0 on U dDm- 

- ■' ’ m—l 


m—l 

Let us assume for the moment that the exterior problem has a unique solution, then we deduce that 

/ M - \ = 

[/ = 0 in K'^\ U Dm ■ Due to the continuity of the single layer potentials, we deduce that C/ = 0 on 

\m=l ) 

M ~ ~ M 

U dDm- In addition, we know that AC/ + k^C/ = 0 for a; G U Dm- From the condition on we deduce 


m—l 


that C/ = 0 in I U D^ 

m—l 


M 


By the jump relations, we have 


m—l 


dU 

— (x) + XmU{x) = 0 
ou 


(KV„)(x) - ^ {V*m, + XmSm,){Oy){x) = 0 


M 


i=i 

jyim 


(3.8) 


and 


dU ~ 

— (x) + XmU{x) = 0 
ou 


(KV™)(x) + ^ {V*m, + XmSmj){Oy){x) = 0 


M 


i=i 

jyim 


for X G dDm and for m = 1,..., M. Here, K* is the adjoint of the double layer operator K, 

f d 

{Kam){x) := / —$k(x, s)cr„^(s)ds, for TO = 1,... ,M. 


(3.9) 


(3.10) 


Difference between (3.8) and (3.9) provides us, dm = 0 for all to. 

M M 

We conclude then that —5I + DL* + DK* + X{L + K) =: —+ P* + AiS : H L'^{dDm) H L'^{9Dm) 

m—l m—l 

is invertible. 


We need now to show that the exterior problem AU + k^U = 0 for x G 


M - 

U D, 

m—l 


with S.R.C 


M 


and + Xm)U{x) = 0 on U dDm, has a unique solution. This result is known under the condition 

^ m—l 


^Observe that adjoint of + T>mm- Hence Fredholm if we show that + T>mm 

is. In [28] , this last property is proved for the case k, = 0 and by a perturbation argument, we can obtain the same results for 
every k. Using the condition on we deduce that — hi isomorphism. 








see 
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for instance. To relax this positivity condition and consider < 0 ior 


M 

5Am > 0 on U dDn 

m—1 

some m’s, we proceed as follows. Set / := (/i,...,/m) with := i(3A)_C/ on dDm where (9A)_ := 
max M __ {—3Am,0}. Hence U satisfies AH + = 0 for cc S Dm \ with S.R.C and + 


u aD„ 

n^l 


TTL — l 


M 


Am + i(?sX)-)U{x) = /m on U dDm- Since now 3(Am + *(5^A)_) > 0, then this last problem has a unique 

m—1 

solution and it can be represented via single layer potentials Taking the normal trace on 

dDm’s, we deduce that (—5I + DL* + DK* + (A + i(3A)_I)(L + K)) where i/i := {ipi, -..'iPm)- But 

from the form of /, we obviously have / = + K)il}. Hence 


(--I + DL* + DK* + \{L + K))iIj = 0. 


(3.11) 


Let us first consider, for simplicity, only one scatterer and assume that the surface impedance is a constant 
A := A*" + i Ah Then (3.111 reduces to 


(_1/ + K*4 

1. If A® > 0, then, as usual in the scattering theory 


= 0. (3.12) 

, we derive that Sip = 0 and then '0 = 0. 
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2. How about A® < 0? Since the operator + K* + A®'S' is invertible in the L'^{dDm) spaces, then the 
equation (3.121 can be reduced to 


(-i/ + K* + X^Sy^SiP = -(zA®)-^ = *(A®)"V 


(3.13) 


i.e. (0,z(A®)“^) is an eigenelement of the operator (—|/ + K* + A’'S)“^S. But (—+ K* + A’’S)“^S is 
compact hence it has only a discrete set of eigenvalues. Then if we take the surface impedance A such that 
z(A®)“^ is different from these discrete values, then ■0 = 0 and hence U = 0. 

We conclude that the scattering by an obstacle with an impedance type boundary condition modeled by 
a constant A®' + zA® is well posed as soon as the imaginary part A® is such that z(A®)“^ is not an eigenvalue 
of the corresponding compact operator (—|/ + K* + X^S)~^S. 


In particular, the operator in (3.12), i.e. —^1 + K* + XS, can be inverted using the Neumann series if A 
satisfies |A®|||(—+ K* + A®'5')“^|| ||f5’f| < 1. This is of course a stronger condition on A but it is enough for 
our purposes. In addition this Neumann series argument applies smoothly to the case where we have variable 
surface impedance’s and multiple scatterers. Indeed, we know that the operator — |l + DL* + DK* + (A + 


M 


z(9A)_I)(L + if) : n L^idDm) 

m—1 

part of A is small so that 


M 

J([ L'^{dDm) is invertible. Hence if the negative part of the imaginary 

m—1 


(5A)_||L + RTIIIK--! + DL* + DK* + (A + z(9A)_I)(L + K))- 


< 1 , 


(3.14) 


then by the Neumann series expansion — + DL* + DK* + X{L + K) is also invertible and hence 0 = 0. 


general condition to inverte (3.11 1 is to assume that 
DK* + (A + 2(^A)_I)(L + X))-1(2(^_(L + K)). 


— 1 is not an eigenvalue of the compact operator {—+ DL* + 









3.2 An appropriate estimate of the densities am, m = 1,..., M 

From the above theorem, we have the following representation of a: 
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a = -(--1 + DL*+DK*+\{L + K))-^{d,, + X)U^^ 


(3.15) 


= -(-^I + DL* + XL)-\I + (-il + DL* + XL)-^{DK* + XK))-^{d,, + X)U^^ 

1 ^ / 1 \ ^ 

= -(--I + DL* + AL)-i ^ ((-2I + DL* + XL)-\DK* + XK) ] (9, + A)C/^", 


if II (—irl + DL* + XL) ^{DK* + AiF)|| < 1. By the assumption a < a %ji/2 1, the operator —^I+ZIL* 

is invertible. We write —|l + DL* + XL = (J + XL{—^I + DL*)~^){—^I + DL*). From the scaling of the 
operators L and (—+ DL*), we show that there exists a constant c, depending only on the Lipschitz 
character of the reference bodies Em’s,, such that if 


|a = max |Am_o|a^ ^ < c, 


(3.16) 


we have ||AL(—+ DL*)~^\\ < 1 and hence the operator —+ DL* + XL is invertible. The condition 
(3.16) is verified if /3 < 1 and a small enough. For the convinience, we denote ||A|| by |A|. 

This implies that 


k < 


|(-il + ZlL* + AL)- 


1- ||(-|I + DL* + AL)-i|| \\DK* + XK)\\ 


\id, + X)U 


In I 


(3.17) 


Here we use the following notations: 

||L>iF* + ATCII := 


{--1 + DL* +XL)-^ 


\DK* + XK\\ , M 

cl n L^idDm), n L-^idDrr 


M 


max 

l<m<M 


max 

l<m<M 


^ c(L'^(dDj),L'^(aDrr^)) 

1=1 

M 


(3.18) 


1=1 

j^m 

{-h + DL* + XL)-^ 


c{ n L^ao^), n ^=(90™) 

= 1 m = l 


max 

l<m<M 


M 

E 

1=1 


{--l + DL* + XmL)-^ 


c(L^(aD^),L^(aD,)) 


max 

l<m<M 


( 2^ Dmm T ^mSmm) 


|cr|l M 

n L^{aD„ 


c(L^(aD,^),L2(aD,r,)) 




:= 

\\U M 

n L^ao^) 

m = l 

= max 1 

l<m<M 

1 \\L^{dD,^) 

and 


:= 

M 

n L^idDr. 

m = l 

= max 

1 ) l<m<M 

1(9 17*11 


(3.19) 

(3.20) 

(3.21) 

(3.22) 


In the following proposition, we provide conditions under which || (—5I + DL* + XL) ^{DK* + AiF)|| < 1 
and then estimate ||cr|| via (3.17). 





















Proposition 3.2. There exists ag depending only on the set of the a priori bounds such that if a < ag aAS 
s <2 — (3, then we have the following estimate 

\Wrn\\i,ll^QD^) < Ce ^ 

where c is a positive constant depending only on the set of the a priori bounds. 

Proof of Proposition 

Suppose 0 < e < 1 and := eB + z C. M". For any functions f,g defined on dD^ and dB respectively, 
we use the notations; 

■= KO ■= + and (g^ix) := g{x) := 9 ■ (3.23) 

Then for each '0 G Lf{dD^)^ we have 

= e^ll'0llL2(aB) (3.24) 

We divide the rest of the proof of Proposition |3.2| into two steps. In the first step, we assume we have a 
single obstacle and then in the second step we deal with the multiple obstacles case. 

3.2.1 The case of a single obstacle 

Let us consider a single obstacle He := eB + 2 with unit outword normal v to its boundary. Then define the 
operator Vd, ■ L'^{dD^) L'^{dDf) by 


(X>z3»(s)=/ (3-25) 

JdD^ Ol'{t) 


Following the arguments in the proof of Proposition 

3.1 

the integral operator —\l + ■ L'^{dD,:) —>• 

Lf{dDf) is invertible. If we consider the problem ( 

1.1 


1 in we obtain 


a={-]^I + Vl^P\SDX\d,P\)U\ 


and then 

\W\\L'^(dD,) < Wi-^I + 'D*d^ + >'SDf)~^\\c(L'^(dD,),L'^(dD,))\\{dv + \)U'‘\\L'^(dD,)- (3.26) 

We have the following lemma, see Lemma 2.4 and Lemma 2.15]. 

Lemma 3.3. Let G L‘^{dDf). Then, 


SD,f^ = e (3.27) 

\\^D^\\c{L'^{dD^),L'^{dD,)) = ^ ll‘^ljll£(L2(3B),L2(aB)) , (3.28) 

(3.29) 

+ + , (3.30) 

+ 4>={{-\l + 'D^^ , (3.31) 











and hence 


--I + V 


D, 


C(L'^iaD,),L^(dD,)) 


—-I + + \Sd, ) V' — 


~2^ + 


--i + v% 


—-I + -\-\eS% ] ij) ] , 


—-/ + + XeS% 


C{L^{dB),L^(dB)) 
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(3.32) 

(3.33) 

(3.34) 


— 2 ^ + 


-1 


C(L'^(dD,),L^(dD,)) 


—-/ + + \eS% 


withS%ip[i) ■■= JgB^''{^,v)'ip{v)dv, 'D'bHO ■= Job and $'^(^,77) := Clc-Jl ■ 


-1 


(3.35) 


CiL^{dB),L^{dB)) 

Ke\^-r]\ 


Let us estimate the norm of 115' 


B\\CiL^{dB),L'^{aB))- 


Lemma 3.4. The operator norm of the compact operator Sd, ■ B{dDf) —>■ L^{dD^), defined in (3.25), is 
estimated by e, i.e. 


.\\C(L2(dD,),L^aD,)) ^ ) ’ 


r(T2(pirt \ r2/'jqn ^ ^ ( 115' 


(3.36) 


Proof of Lemma 3-4 To estimate the operator norm of Sd,, we decompose So, =■ 5)5 = S]^ + Sff into 
two parts ( independent of k ) and ( dependent of k ) given by 


SdMx) ■= 


1 


iao, 47r|a; - y\ 


'fp{y)dy, 


r iK\x-y\ _ 1 


(3.37) 

(3.38) 


With this definition, 5|^ : L'^{dDf) —>• L'^{dD^) and 5^" : L'^{dD^) —>• Lf{dD^) are compact. From (3.28), 


it can be observed that 

Il5' 


_ , II C'i’K. II = f II II 

DA\c(L'^(aD,),L^(dD,)) ^ll‘-^B \\c(L‘^(dB),L‘^{aB)) ^ B\\c{L'^{aB),L'^{aB)) 


(3.39) 


On the other hand, as mentioned in (2.30) of 15, Lemma 2.5], the following estimate can be obtained, 

1 


5 


D, 


< —Ke^|i95| for Kynaxdiam(D^) < 1. 

c{L^{aD,),L^{aD,)) 2 t: 


(3.40) 


Hence the result follows. 


3.2.2 The multiple obstacle case 
A way of counting the small scatterers 

Before proceeding further we make the following observation. For m = 1,..., M fixed, we distinguish between 
the obstacles Dj, j ^ m by keeping them into different layers based on their distance from Dm- Let Tim, 






























1 < m < M he the cubes of center Zm such that each side is of size (| + c?“) with 0 < a < 1 and it contaiJi^ 
only Dm- Let us suppose that these cubes are arranged in a cuboid, for example unit rubies cube, see FigUl 
in different layers such that the total cubes upto the layer consists (2n + 1)^ cubes for n = 0 ,..., [d““], 
and ^Im is located on the center. Hence the number of obstacles located in the n ^ 0 layer will be 
[{2n + 1)^ — (2n — 1)^] and their distance from Dm is more than nd°‘. 



Figure 1: Rubik’s cube consisting of two layers 


Lemma 3.5. Form,j = 1,2, the operator Smj ■ L^{dDj) —>■ L^{dDm) defined in Proposition 

see (3.5), satisfies the following estimates, 


3.1 


For j = m, 


\\^mm\\c{L'^(dDrr,),L'^(dDrr,)) — ^ 


5^ 


+ ) I 


(3.41) 


for Kme^xO < 1. 

For j m, such that Dj G Nff, n = 1,..., [d““] 


15, 


< 


1 1 


™Jll£(L2(9Z>,),L2(az3„)) ^ 47rnd“ 


m\' 


(3.42) 


where \d]d\ '■= max|9i?m| ond 


S 




__ 11 11 


Proof of Lemma 3.5 The estimate (3.41) is nothing else but (3.36) of Lemma 3.4 replacing B by Bm, 
z by Zm and by Dm respectively. The proof of the estimate (3.42) is a straightforward consequences of 

□ 


(2.37) in 15 Lemma 2.6]. 


Proposition 3.6. Form,j = 1,2,...,M, the operator : Lfi{dDj) —)• L‘^{dDm) defined in Proposition 


3.1, see (3.6), satisfies the following estimates, 


For j = TO, 




-1 


5 C^mj 


(3.43) 


CiL^(dD„,),L^{aD„,)) 


where C^m '-= 






27r-K2e2|9B^| 


\\c(L'^(dB„,),L'^{dB„,)) 






























For j ^ m, such that Dj € 7V^, n = 1,... ,[d “] 
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< 


1 / K 


mj\\CiLHdD,),L^dD^)) - 4 ,^ n^d^' 


1 


15^1 


where \dl^\ := maxSBm- 


In addition, as a consequence of (3.34), we can also prove that 
• For j = m, 


(3.44) 


where C^m '■= 


2 ^ + T^mm F )^mSmm 






C{L^{dD„,),L^{dD,^)) 


I £(1,2 (8B„,), 1,2(83^)) 


27r-(|A,„|+K)«;e2|aB H 

'll'- 2 Bra Bra/ 11 £ (1,2 ) , 1.2 (aB„ )) 


(3.45) 


Proof of Proposition 3.6 This result can be proved in the similar lines of the proof of 15 Proposition 

2.17]. 

□ 


End of the proof of Proposition [gTl . By substituting (3.42) and (3.44) in (3.18), (3.45) in (3.19), 
and using our discusion related to Fig on how we count the number of small scatterers, we obtain 


and 


M 


\DK* + \K\\ = max ^ 

l<m<M 


™Jlk('C2 


i=i 

j^m 


C(L'^(dDj),L^(dDa,)) 


[d- 


< Y.[{2 n+lf-{2n-lf\- 

n—1 

M-“] , 


[d 


1 ^ |A| + K ^ 1 


47r V u(i“ n'^d'^° 


|9B|, 


^ ^[12n2 + l] 


|A| + K 

1 1 ^ 

1 |9B| e' 

nd°^ 

n'ifp.a ) 

'|A| + k 

1 ' 

)|5B|e^ 

nd°^ 

n^d'^'^, 


1 

27r 

1 


[d I , Id “] 

(|A| + K)d “ ^2 [12u H ] + d ^ [12 H -] 


n—1 


n—1 


|5B|e^ 


< — [(|A| + K){6d-^°‘ + 7d-2“) + 13d-3“] |aB| 


(3.46) 


--I + DL* + AL 


-1 


= max 
l<m<M 


2^ '^rnm ^m^mm 


-1 


C{L^{^Drr^),L^{^D^)) 


= max Cem- 
l<m<M 


(3.47) 


Hence, ( 3.46|3.47 ) provides 


(--I + HL* + AL)-i IIZHF* + AiF|| 








































< ( max Cem ) |i9B| ^ [(|A| + k )(6 + 7fi“) + 13] d . (3.48) 

\l<m<M / zTT 

= :C, 


By imposing the condition ||(— |l + HL* + XL) ^|| \\DK* + XK\\ < 1, we have from (3.17) and ( 3.20|3.22 ); 
Il^"i|li2(a£)^) < ||o 


< 


I|(-il + iAL* + AL)-1|| 

||(5, + A)17^" 

1 - 

||(-il + OL* + AL)-i 

IIHAT* + aa:|| 

< 

Cp 

(-il + 74L* + AL)-i 

max (9,. + A)17* L 2 -„„ . 
l<m<M "' (on^) 

< 

C max \\{di, + X)U'^\\ ^. . (c\=Cp 

l<m<M " {dO-ra) \ ^ 

max Cem ) 
l<m<Af J 


(3.49) 


for all TO S {1, 2,..., M}. But, 

||(9j^ + X)U ||j;,2(a£)^) < ||t^ Wl^od^) “*" \\L^(dD^) 

= |A|e +fce|9B„|5 (Since , C/*(x, 0) = 

< (|A|+K)e|5B|5 ,Vto = 1,2,...,M. 


Now by substituting (3.50) in (3.49), for each to = 1,..., M, we obtain 

ll^m|lL 2 (a£i^) < C{k)€, 

where C(k) := C |9B| = (| A| + k). 

The condition \\{—^l + DL* + AL)“^|| \\DK* + A7ir)|| < 1 is satisfied if 

Cs = ( max Cem^ |9B| [(IA| + k)( 6 + 7d“) + 13] cr^“e^ < 1. 

V l<m<M ) ZTT 


(3.50) 


(3.51) 


(3.52) 


Since Am = AmO « ^ and d « o‘, in particular d > dminCi*, then (3.52) reads as a 3“*+^ ^ < c, where we set 

-1 


C •— ( [('^+ + 'tmaxO^)(6 + 7[finiax] ) + 13a^] 


1 


|5B| 


27r [ max diam(Bm)]'^ i<m<M 
l<m<M 


max 


> 1 (3.53) 


with A+ := max lAmol and with the rewritten form of Cem mentioned in Proposition 

l<m<M 


3.6 


as 


27r 


Cem ■— 


_ir_l_70*"* J_ -^mpe ^ ciK. 

2^^ ' ^Bjri ' f max diamiBm)]^ 

l<m<M 


C{L^{dBm),L^{dBm)) 


'■l<m<M •' 


_ir-LT)*'«* J_ ^ 

2 Bm [ max diam(Bm)]^ Bjj 


-1 


Observe that s = 3at. Hence (3.52) makes sense if s < 2 — ,5 and A+ satisfies (3.53). 
Again, since Am = XmooT^, (3.51) can be rewritten as 


I11^2(31)^) < C{k) 


-i-t 


C{L‘^{dB^),L‘^{dB^)) 


(3.54) 
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C{k) 


c|aB|M- 


max diamiBr, 
l<m<M 




|^<C|9B|5( 


[ max diam(Bm)\^ 
l<m<M 



□ 


3.3 Approximation of the far-fields. I. Approximation by the total charges 

We start with the definition of the total charges Qm: m = 1,... M. 


Definition 3.7. We call the (Jm’s used in (3.1), the solution of the problem (1.1 1.3), the surface charge 
distributions. Using these surface charge distributions, we define the total charge on each surface dDm 
denoted by Qm as 


Qm • 


IdDr, 


i{s)ds. 


(3.55) 


In the following proposition, we provide an approximate of the far-fields in terms of the total charges 


Qm,' 


Proposition 3.8. The far-field pattern U°° of the scattered solution of the problem ML has the following 
asymptotic expansion 

M 

u^{x) = + 0^-/3)]^ (3.56) 


with Qm given by (3.55), if a < 1 where 0{Ka^ < Cko? ^ and C := . 


max diam(Bm) 

l<m<M 


Proof of Proposition 3.8 From (3.1), we have 


M n / \ 

U^{x) = y^ ^,,{x,s)am{s)ds, for X eR^\( U Dm] ■ 


Hence 


M 


U^{x) = Y. 

m—1 

M 

= E 

m=l 

M 

= E 


IdDr, 




e-i-'S:--rrr^mis)ds+ / [e- 


Wmis)ds 


'dD„ 


IdDr, 


-IKX- Zrr 


Qr. 


IdDr, 


^]amis)ds . 


(3.57) 


It is imp ortant to remark that, if we do not dist inguish the near by and far obstacles, as it is discussed in the beginning of 


section 


3.2.2 


and by following the way it was done in 


15 


16 


t get the estimate \\DK* + A/C|| < ^ 1^91 


in place of ( |3.46[ l and hence the condition | |3.52| l will be replaced by (M — 1) < ^0? for some suitable constant cq. However, 

this condition is too strong to enable us to apply our asymptotic expansion to the effective medium theory where we need to 
choose M ~ with s = 2 — /3 and d ~ a^with f > f. 
























For every to = 1, 2,, M, we have from Proposition 3.2 

\(Jm{s)\ds 
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>dD„ 


< l|l||L2(aD„)|km||L2(g£i^) 

I|1||l 2(9D„)C |5B|= (A+ + Kmax) 6^“^ 

< \dB\C{- 


, 2 -/: 


with C :== 


max diam(Bm)]^ 

"l<m<M 

|9B|C(A+ + Kmax[ max diam{B,n)]^) 

l<m<M 


max diamiBr^ 

l<m<M 


It gives us the following estimate; 


ldD„ 


^ ^’-]am{s)ds 


< 

/ 

JdD^ 


< 

j 

JdD^ 

oo 

tvls — Zr 

1^1 

< 

J 

JdD^ 

oo , 

i:»'(f) 

Z=1 

< 

(|3.58t 

Ca^- 

Z=1 

_ 

^Ck 



i{s)\ds 


“ K 


which means 


r fora< 

JdDrr, 


(3.58) 


(3.59) 


(3.60) 


(3.61) 


Now substitution of (3.61) in (3.57) gives the required result (3.56). 


3.4 Approximation of the far-fields. II. Estimates of the total charges 

3.4.1 Derivation of the linear algebraic system 

We start with following a priori estimate on the total charges Qm, to = 1,..., M. 

Lemma 3.9. For to = 1,2,..., M, we have 

\Qm\ < (3.62) 


where c := \dld\G{- 


max diamiB^ 

l<m<M 




Kinax) with did and C are defined in (3.42) and (3.49) respectively. 


Proof of Lemma 


3.9 


The proof follows as below; 



<7m{s)ds 


\Qm\ 





















< l|l||L2(aD„)lkm||L2(9D„) 

< IMII 
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I|1||l2(9d„)C |5B|2 (■ 


< |(9B|C(- 


max diamlBrn)]^ 

l<m<M 

, ^ .2-/3 


I ^ W1-/3 

~ '^max / 


max diamiBrri)]^ 

l<m<M 


□ 

The following proposition gives an approximate characterization of the total charges Qm^ m = 1, 

Proposition 3.10. For m = 1, 2,..., M, the total charge Qm on each surface dDm of the small scatterer 
Dm can he calculated from the algebraic system 


M 


0^ = -U\Zm)-^Cj^^{Zm,Z,)^+ErT 

Go 


a 


i=i 

j^m 


(3.63) 


where Err := O ^ and Cm ■= —^m\dDm\- 

Proof of Proposition 3.10 For Sm € dDm, using the impedance boundary condition (O , we have 


n _ !„ \ \ \ TTtr„ ^ _ ^m{Sm) 

0 ('^m) T ^mU (^TTi) — 

avm. 2 


IdDrr, 


M „ 

{sm,s)am{s)ds + {sm, s)aj{s)ds 

JoD. OVm 


, = 1 JOD, 
j^m 


^ ^ r gjji 

“t ”I (^m; ^^O'mi^^ds p — ('^m) “t ^raU (^m) 

JdD, dVm 


i=i 

Integrating the above on dDm, we can write it as 


(3.64) 


1 [ CrmiSm)dSm+ [ ( f {Sm, s)dSm] Crm{s)ds + y^ [ ( f {Sm, s)dSm] CTj {s)ds 

2 JdD^ JdDrr, \JdD^ dVm J ^ Jan, \JdD^ dvm J 


j^m 


M 


+Am yy^ 


j = l-'dD, \JdD^ 


^Ki^m,^)dSm I ^m{^)ds — 


r dw 

I dDm 


{^m)dSm j ^mC\sm)dSj., 


'dD„ 


It can be rewritten as 

f ( f d% 

IdDm \JdDm 


2^" 


(Stti, s)(iSj7j j Cm {s)ds+yy [ ([ ^{Sm,s)dSm] crj{s)ds 

/ _, j^m'^dD.KJdDmdVm J 


= :A 


=:B 


+ Am / ( / 4>„(sm,s)dsm ) crm(s)(is+Am ^ / ( / ^fc(sm, Zj)dsm j crj(s)ds 

JODm \JdDm J \JdDm J 


= :C 


=:D 


r dw 

I dDm 


(Sm)dSm I ^niC {Sm'jdSm P A. p XmD , 


I0D„ 


with 


A' := 


I dDm \ddD 


dl^m dl^m 


dsm ) crm{s)ds 


(3.65) 

























M 


D .— ^ ^ '5) ^^j')] ) CFj{^s)ds. 


We can approximate A and A' as follows; 

^ = [ 


d^o 

dvr, 


(Sttx, s)(^5j7j j (7-fYi {s)ds 


IdDr^ \JdD^ 

lK^£^{l)]{s)am{s)ds 


ldD„ 

1 


^ JdD^ 
^Qm- 


arn{.s)ds 
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(3.66) 


(3.67) 


Here is the double layer operator defined as in (3.10) but with zero frequency and on the boundary 
of Dm- Observe that, 


|H'| = 


< 


'dD,n \JdDrr 

f (/ 

IdD^ \JdD„ 

/ (/ 

IdD^ \JdD„ 

O . 


s) _ (g s) 

UUm 


47r\Sm - sP 

9 oo / \ 

1=0 


E( 

1=1 


dSm ) (Jra{s)ds 

1 


lK\Sr„. - S\ 


\l 


ll {1 + 1)1 


dSm |crm(s)|ds 


Now, we can approximate B 


M 

^ = E 

j^m ^ 
M 

= E 

j^m ^ 
M 

= E 

j^m ^ 
M 

= E 

[d-“] 


{Sjn,s)dSm ) crj{s)ds 


dDj \JdDrr, ^’^rn 

[ (^J A$«(?/ m 5 s)d2/m^ o'j{s)ds 

[-K^^f,{ym,s)]dyjn ] crj{s)ds 


'dDi \JD. 


IdDi 


O K 


aj{s)ds 


^[(2„ + 1)3_(2„_i)3] 

5-P [d”“] 


n—1 


o 




nd‘^ 


= 0\2k^ 


= O 2k 


d° 


E[i2^ 

n—1 


l+d + 7] ■ 


dsjYi j C7jyi{s)ds 


(3.68) 


(3.69) 
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Since s) — Zj) = O {-^ + for s € dDj,j ^ m, we can approximate D and D' as 

follows; 


^ = E 

j/m 

M 

= E 


^KiSm,Zj)dSm (7j{s)ds 


'dDi \JdD 


dDi \JdDr, 


^:^{z„i,Zj)dSm <Jj{s)ds 


IdDj \JdD 


Zj^ ^ZjY^ds 

m I {s)ds 


= ^ ^K{zm,Zj)Qj\dD^\+0 (^K+\dDm\ J aj{s)ds'^ 
“ r ^a^~^ ( 1 w 

= 'y ' ^K{Zrm Zj)Qj \dDm \ + O I —- ( K + —- | | 

L \ “mj V Umj J J _ 

M /n^-P ( 1 

= ^ ^f^{zm, Zj)Qj\dDm\ + ^ [{2n + 1)^ — (2n — 1)^]0 f f k + 

j^m n—1 ^ ^ 


^ I n^-P ^2■n'^ -I- 1 / 1 

= ^f^iZm, Zj)Qj\dDm\ + O { 2 —^ ^ 

y n=l ^ 

Jl£, / o5-/3 r fiK+ltlX 

= + O ^2 7kH J 


li’i = H 




('^mi'S) -^j)] (Tj(5)ds 


“ E^(^('"+^ 

\ ^rnj \ ^mj 


n—1 ^ ^ ' ' 


L 6s;+ 13 


= O 2 -— 7k+ — 
V d?°‘ d' 


Now, let us approximate C. Since |$K(sm,s)dsm| < 4^ Jqjj^ = 0{a), then 


\C\= [ ^KiSm,s)dSm] Crjn{s)ds 

JdDrr, KJdDrr^ / 


= 0 {a^-^) 


(3.72) 



Hence, from ( 3.65|3.72 ), we obtain the approximation below 
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M 


Qm H" ^ ^ ^K{Zrm Zj)X m \dD 

m \Q, 

j^m 

= + O ((|A™| + K)Ka^) + X^O 


+0 {k^g^ + 0 i + 7] ) + X^O ( 2 


.5-e 


d2c 


7k + 


6k + 13 


(3.73) 


Indeed, Xm{{U''{sm) - U^{zm))dsm = 0(|A|Ka^) and §^{sm)dsm = 0{K'^a^). In addition, since 
/3 < 1, then A^O {a^-d) = 0{a^-'^d) = o{a'^-d) = o(A|5£i^|e*''’®-^™)‘. 

We can rewrite the above algebraic system as 


1 

Xm\dD„ 


M 


.Qm = -e*''®'*’" - ^ ^K.{.Zm,Zj)Qj + Err 


j^m 


with 


Err := 


1 


Xm\dDr, 


+0 2 


75 - 2/3 


d2c 


(('^++ ^) + O (k^o'^ ^) + O ^2k^ ^2q, + 7] 

+ O (a3-2/5) 


7k + 


6k + 13 


I 


A,„|aL»„ 


O a^-^d + 


^5-2/3 


= O + 


.3-/: 




In the last two lines in the above approximation, we used the fact that k < Kmax and d < d„ 


(3.74) 


(3.75) 


3.4.2 Invertibility of the algebraic system 

We define the following algebraic system 


M 


Or ' 


j/m 


c. 


for all TO = 1, 2,..., M. It can be written in a compact form as 

BQ = O, 

where Q, G and B G are defined as 


( 


1 

Cl 


-$k(zi,Z2) -$k(zi,Z3) 


B := 


-^ k ,{ z 2 , Zi ) -$^(22,2:3) 

V -$„(zm,-Zi) -$k(2:m,Z2) ■■• -$k(zm,-ZM-i) 

0 := ( <9i Q2 Qm and O := ( C/*(zi) C/*(z2) 


-$„(zi,zm) \ 

-$k(2 ;2,2 :m) 


/ 


U^{zm) 


(3.76) 


(3.77) 


(3.78) 


(3.79) 


The above linear algebraic system is solvable for Qj, 1 < j < M, when the matrix B is invertible. Now we 
give sufficient conditions for the invertibility of system (3.77). We recall that Am := Xm,oa~d, /3 > 0. 























Lemma 3.11. We distinguish the following two cases: 
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Let o) < 0 r o-nd assume that min then the matrix B is invertible and the 

y l<m<M 


TT dt 


solution vector Q of (3.771 satisfies the estimate 


M 


E l0mP<4 


Let 5ft(Am,o) > 0 and assume that 


min 3?C„ 

l<m<M 


-2 


V2M^, 

max ICmP TTtit 

1<m<M 


M 


E ■ 


(3.80) 


min 3ff(—Cm.) _ 

r^~ - rp^, —TT^ > ^then the matrix B is invertible and the 


solution vector Q of (3.77) satisfies the estimate 

min 3?(-(7m) 


M 


E ic 


l<m<M 




M 


max ICmP 7 rdt 

l<m<M 


E \U\^m)\ ■ 


(3.81) 


where M^ax '■= M a®, recalling that M := 0(a ®), as a ^ 0. 
Since Cm '■= —^m\dDm\ and d > dminof, the condition 


min |5RC„| ._ 

jg satisfied if 

TT dt 


max ICrnP 

l<m<AL 


° ,2 -- > This is possible if s < 2 — 13 and a is small enough or s = 2 — /? and A_ and 

ndt 

A_|_ satisfy Recall that t > hence f < 3. If we assume dmin C 1, then the last condition 




is satisfied if and if dmin > 1 then we take ^ proof of this lemma will be 






given in the appendix. 


We can rewrite the inequatilies (3.80) and (3.81) using norm inequalities as 

-1 

VMl, 


M 


E ic 


< 2 


m—1 

M 


min iRC„ 

l<m<M 

max I Cm 

l<m<M 


max I C„ 


7T l<m<M 


M max \Cm\ max \U'^(zm)\ (3.82) 


E 1C 


min 3?(-C'm) 

l<m<M V^-^? 


-1 


m—1 


max \Cn 

l<m<M 


, max \Cm\ M max \Cm\ max {7®( 2 ^) (3.83) 

Trdt l<m<M ; l<m<M l<m<M ' ' " I' 


which holds for the cases 3?Am < 0 and 3?Am > 0 respectively. 


3.4.3 The dominant part of the total charges 


The difference between (3.76) and (3.63) produce the following 

M 


Qm Qvi 


= - E (C - Qj) + Err, 


(3.84) 


1=1 

j^m 


for TO = 1,... ,M. Comparing the above system of equations (|3.84|) with (|3.76) and by making use of the 

(3.85) 


estimates (3.82) and (3.83), we obtain 

M 

E “ Q™) = O {Ma‘^~^Err) . 


m—1 


In this case, we can actually relax the condition s < 2 — /B on the number of small scatterers s, see Remark 6.1 













































We can evaluate the Qm’s from the algebraic system (3.761. This means that Qm’s are the dominant pai^l 
of the total charges Qm’s- 


3.5 Approximation of the far-fields. III. End of the proof of Theorem |1.2 

Using ( |3^ in ( |3.56[ ) we can represent the far-field pattern in terms of Qm as follows 


M 


m—1 

M 

= Y. [Qm + (Qm -Qm)]+ 0 (« a^-^)] 

m—1 

M 

= ^ e-*""-"™0^ + O(Ma2-'5[Urr+(Ka)]) 

m—1 

^ \ 

^ + O ^ j 

M 


(3.86) 


since, as s = 3ta < 2 — /3, we have ^ ^ = 0{a) = o(a^ if /? < 1. Hence Theorem 

proved with the replacement of in the statement by C„i- 


1.2 


IS 


4 Proof of Corollary |1.3 


The steps of the proof of Corollary |1.3| are the same as for the proof of Theorem Here we only explain 
the main changes that are needed to derive it. 

We recall that, for an obstacle of radius e, S{(j>){s) := and 'D{(j)){s) := 

IdD, Similarly, we set 5 g((/>)(s) := G^{s,t)(j){t)dt and T>g(</')(s) := ^%j^(l){t)dt. 

We see that lU„(a;, z) := Gk,{x, z) — z) satishes 


(A + k^u^)1Uk = K^(l — in ( 


(4.1) 


with the Sommerfeld radiation conditions. Since $k(-, z), z is bounded in LP(fl), for p < 3, by interior 
estimates, we deduce that W{-,z), z G is bounded in for p < 3, and hence, in particular, the 

normal traces are bounded in L?{dDf). Then we can show that the norms of the operators 


and 


Sg-S: L^idD,) H\dD,) 


Vg-V: L'^idD,) -)■ L^{dD^) 


are of the order 0(e) at least. The representation (|3.1|, in Proposition 3.1 needs to be replaced by 


M 

Ui{x) = V*ix) +Y G^{x,s)am{s)d. 
m=l -^90^ 


.,, , M - 

S, X G M'^\ U Dm , , 
m—1 


(4.2) 

(4.3) 

(4.4) 


where is the total field corresponding to the background modeled by the index of refraction n, see (1.151. 


We use the single layer potentials defined by the Greens’ function instead of the fundamental function 



















The main tools used in justifying Proposition 3.1 are the invertibility properties of the correspondiSf 
integral operators (i.e. the Fredholm property) and the jumps of the double layer potentials defined by <1 ?k- 
These two tools are satisfied also when we use instead of due to the error bounds in (|4.2|) and (|4.3). 
The a priori estimates on the densities a„ 


I’s derived in section are quantitative versions of the result 
in Proposition 3.1 Due to error bounds in (4.2) and (4.3), those estimates can then be translated to the 
densities used in the representation (4.4). Of course, in Lemma 3.3 one needs to replace the equalities by 
inequalities to estimate the properties of the operators defined by on the scaled obstacles D^, in terms 
of the properties of the operators defined by on the original obstacles B. 

Finally, to do the same analysis as in section 3.3 one needs to split G„ as G^ = $„ + (G^ — $„) and use 
the results derived in section 3.3 and again error bounds in (4.2) and (4.3). 


5 Justification of Remark 1.4 


We show only the main changes needed in the proof of Theorem m This change occurs in the proof of 
Proposition 3.10 and precisely in evaluating the term G, i.e. (3.72). We rewrite G as 


G — / ( / J (T'^i^S^ds F / I / s) ^')dSm J (T^(s)ds. 

JaDrr, yJdDrr^ J ddD^ \ddD^ J 


Since 


/ ( / '5) j (JjYii^S^ds 

JdD^ \JdD^ J 


1 ^ Ac'a'-i 

Z=1 


\dDm\ / |crm(s)|ds 

JdD^ 


= O {Ka‘^-1^) 

then we can write G as follows; 

G = / [ / $o(sm,s)dsm ) crm(s)(is + O 

JdDrr^ \JdDm J 

= ImQm + O , 

where Im = ^o(sm,t)dsm, is a constant for each t G dDm if is a ball. Indeed, for the sin¬ 

gle layer potential >) (s) = JoD, <i>o{s,t)^it)dt, we have the jump condition as (^S]^ t/j) (s)|_ = 


(5.1) 


I + K 


tf}{s) recalling that is the adjoint of the double layer operator. Since [1] = [1] = 

— 1/2, then we have [1] = 0 on dDm- Hence [1] is a constant, namely Im, as it satisfies 

A [1] = 0 in Dm and has zero normal derivative on dDm- 


With this correction at hand the estimate (3.73) becomes 


M 


( If ^ml-m^Qm F E m \dD 

m \Q, 

= F O ((|A^| F K)Ka^) F A^O (a^"^) , 


5-t 


+0{K^a‘^ +0 [2 k‘^-^[—+ 7]j + XmO (2 — 


and the system (3.74) as 

1 ^m^rr 


M 


Xm\dDr, 


■Qm = ^K{zm,Zj)Qj +Err 

j=/=m 


7k- 


6k F 13 


(5.2) 


(5.3) 































with 
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Err := 


1 


\m\dDr, 


0 ((A+ + «»)»»-«) + O + O + 'l) 


+0 2 


, 5 - 2/3 


<P^ 


Ik + 


6k + 13 


+ O (a4-2/3) 


A,„|9D„ 


O a'^-P + 


j5-2/3 

d3“ 


— O I a + 


. 3 -/; 


d3° 


(5.4) 


Finally, the estimate (3.861 becomes 


U°°{x) = 


M 

^[e-*«"'^’"g„ + 0(Ka3-/5)] 

m—1 

M 


= Y. [Q™ + {Qnr - Qm)] + 0 {k 

m—1 

M 

= ^ e-*""'"'"Q„ + 0(Ma2-^[£;rr+(Ko)]) 

m—1 
M 


(O) 


-IKX-Zrr 


m—1 

M 


Qm + O I Ma I a + 


a 


3-t 






m—1 


(5.5) 


3-/3 


since, as s = 3ta < 2 — we have ^ 3 ^ = 0{a^ ^ = 0{a). 


6 Appendix: Proof of Lemma 3.11 


We start by factorizing B as B = —(C~^ + B„) where C := Diag{Ci,C 2 , ■ ■ ■ ,Cm) € I is the 

identity matrix and B„ := —C“^ — B. We have B : —)■ so it is enough to prove the injectivity in 

order to prove its invetibility. For this purpose, let X, Y are vectors in and consider the system 


(c-i + B„)x = y. 


( 6 . 1 ) 


Let (.)’'®“* and (•) 
ber/vector / matrix. 

(Ip); 


denotes the real and the imaginary parts of the corresponding complex num- 
For convenience, let us denote C“^ by C/. Now, the following can be written from 


(Cjea/ grea/^^reai _ ^(.^mg ^ ^img^^^img 
(C^ea/ ^f^^mg ^ ^img^^^real 


Y 


real 


Y'^mg 


( 6 . 2 ) 

(6.3) 


which leads to 


((c7“' + b;“')V'"9,x™ 3)+((c7® + B™s)x’'“',V™9) = ( 6 . 5 ) 

















By summing up (6.4) and (6.5) will give 
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' ^real\^ _j_ ^ grea/^rea/ ^real\^ _j_ ^ Q^real^img _j_ ^ -Qreal j^img j^img'^ 

_ ^ yrreal j^real'j _|_ ^ yrimg j^img'j 


We can observe that, the right-hand side in (6.6) does not exceed 


^j^real j^real'^l/2 ^Yreal yrreal'^l/2 ^ j^img j^img'^1/2 ^Yimg -Y^img'sXj 2 

At this stage, we divide the proof into two cases. 

1. If 3?(Am,o) < 0 for each m. In this case 3?(7m > 0. We know that: 


real ii i \rreal \2 




where := and := if * J and 


i,j = I, ...,M. Hence | i ^ j- Arguing as in (|3.46|), we see that 


M 


-2a [1^1 2M(i-3“ 


E (sr'jfj < *f E 1(2" +1)’ - (2" - < «— E1 = 


j,i=i 


Observing that M = 0{a ®) < MmaxO, ® and that s = 3ta, we obtain: 




From (6.6) and (6.7), we deduce that 


min 3?C'm , 

l<m<M \/2ilAr 


M 


M 


( max ICml)^ Trdt 


^ < 2 ^ |A„ 


1/2 / M \ 1/2 

I ( X! J 


V,m=l 


yrn=l 


which yields 


M 


E < 4 


min 3fi(Cm) 

l<m<M y/ZMyy 


M 


m—1 


( max \Cm\Y Trdt 

l<m<M 


E 


m—1 


Thus, if 


min Sf(C„) 

l<m<M 

( max ICttiD" 
l<m<M 


> 


V2Ju 


r d t 


then the matrix B in algebraic system (3.77) is invertible. 


2. If 3fi(Am,o) > 0 for each m. In this case 3?Cm < 0 and so to prove the invertibility of B in 
system (3.77), consider (—C“i — B„)X = E in place of (6.1) and proceed in the way as it 
the case SliAm > 0 for each m. Then we can get the following estimate 


M 


E <4 


min 3fi(-C'„i) 

l<m<M 

( max \Cm\)^ 

l<m<M 


^/2Mjnax 
7T di 


-2 

M 


Ei^™i'- 


and the invertibility of the matrix B, under the assumption that 


min m-Crr 
l<m<M 

(max I Cm. I) 

l<m<M 


> 


V2m:„ 


( 6 . 6 ) 


(6.7) 


( 6 . 8 ) 

:= 0 for 


(6.9) 


( 6 . 10 ) 


( 6 . 11 ) 


( 6 . 12 ) 


algebraic 
done for 


(6.13) 
































3. We are left with the case where 3?(Am,o) > 0 for few m’s. In this case we multiply every line of t8.§ 
system (6.1) corresponding to 3?(Am,o) > 0 by —1. Hence, the system (6.1) becomes 


(C^i+B„)X = Y. 


(6.14) 


where now every component of the diagonal matrix C ^ is positive and ||B „||2 = ||B„|| 2 . Then we are 
in the case (1). 


□ 


Remark 6.1. Assume that 5i(Am,o) < 0- Following the computations in 15 Lemma 2.22] and assuming 

min 9?Cm 


that 


571 l<m<M _ 

3 ( max \Cm\)'^ 

l<m<M 


> ^ and 7 := min cos(K|zm — Zj\) > 0, we can prove that 




min iRCm 

l<m<M 

( max |(7^|)2 

l<m<M 


37 
Stt d 


M 


M 


E < 2 E 


m—l 


\m—l 



which yields 


M 


E ^ 4 


min iRCm 

l<m<M 

( max \Cm\)‘^ 

l<m<M 



-2 

M 

E 14^-1'- 

771=1 


and thus the invertibility of the matrix B in the algebraic system (3.77). 

min UiCm) , _ 

jg satisfied ift<2 — fi. 


( max I Cl 

l<m<M 


. 1 )" 


Observe that the condition 


We can see that if the number of the scatterers is 0{a~‘) such that s <2 — (3, then we do not need upper 
hound on t, which allow as to have very close scatterers. However, if we want to have larger number of 
scatterers, i.e. s limited only by the hound s < 3, then we need a condition on t, i.e. t <2 — fi, which makes 
more restrictions on the minimum distance between the scatterers. 
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